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Abstract 

Rigid contact model for course-grained soil particles is 
established combined with the distribution of contact force 
and normal of contact force based on the research of contact 
characteristics of particles and ignoring the deformation of 
particle. The local constitutive model is acquired after 
analyzing the relationship between contact force of particles' 
and local stress in course-grained soil. Furthermore, a two 
dimensional discrete constitutive model is built up. It is 
proved through discrete constitutive model that the fabric 
will change during the deformation of course-grained soil 
which results the change of physical and mechanical 
characteristics. The change of fabric affects the macro 
mechanical responses characteristics of course-grained soil. 
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Introduction 

Coarse-grained soil such as gravel soil, soil-rock 
aggregate mixture and rockfill is a typical granular 
media composed of discrete particles within two or 
three orders of magnitude in dimension. And it is 
widely applied in highway, railway, hydraulic dams 
and high embankment engineering. The physical and 
mechanical characteristics are quite different from 
traditional solid medium. Coarse-grained soil is a 
loose media in substance, and the contact between 
particles is discontinuous point-contact. From the 
paper of Luan Maotian, for coarse-grained soils and 
other loose granular media, the research of physical 
constitutive model of deformation and strength 
characteristics to overcome inherent limitation of 
mechanical methods of continuous medium with basis 
of micro mechanics is quite potential. Therefore, 
particle characteristic of coarse-grained soils, 
interaction and essential fabric of spatial uniform 
distribution between particles can be descripted 


objectively. Based on mathematical function of 
statistical distribution of the reflected fabric 
parameters, the relationship could be set up between 
inherent evolution of microstructure fabric parameters 
coarse-grained soils and macroscopic mechanics 
responses . 

Granular media mechanics is assumed that coarse- 
grained soil is composed of solid particles contacted 
with each other, and interaction between particles 
obeys the laws of probability. Granular media 
mechanics is applied in studying the mechanical 
phenomena on the contact points of particles, and 
describes the phenomena according to the formula of 
mathematical statistics. Some mechanical models have 
been built by deferent research. However, those 
models can't interpret the relationship between fabric 
change and mechanical response of coarse-grained 
soils under load very well. 

Rigid Particle Contact Model of Coarse- 
Grained Soil 

There is a certain gradation within coarse-grained soil 
particles where large particles act as the skeleton and 
small particles fill in the void between the large 
particles. Both the interaction between particles and 
the filled surrounding medium and particles are 
belong to solid contact mechanics in terms of the 
mechanical characteristics. It is supposed that the 
number of coarse-grained soil particles is enormous 
amount, so the macro mechanical parameters have 
statistical significance. The uneven stress of 
microcosmic can be described by average stress. When 
the coarse-grained soil particle number is infinite, and 
coarse-grained soil is continuous in the certain space,, 
the summation in the calculation could be changed 
into integral. 

The particle contact normal density distribution 
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function E(n) is introduced, where n is the contact 
unit normal vector. The number of contact points 
within the n—>n+dn is E(n)dn . Supposing f(x a ,n) 
is the contact force component located at the contact 
point x a with normal n . For a certain volume V of the 
coarse-grained soil, the total contact vectors is zero 
based on the conditions of static equilibrium. 

\fi(n)E(n)dn = 0 (i = x,y) (1) 

Where f,(n) is the average of contact force within 
n + dn . 

The contact between particles is one of the most 
fundamental problems in particle mechanics model. 
Assumed that the coarse-grained soil particles are 
rigid, and the deformation of particles is ignored, so 
the contact of coarse-grained soil particles is elastic 
contact. Coarse-grained soil appears sliding or rolling 
possibly under load, so that two adjacent particles can 
be viewed as two rigid bodies connected by deformed 
springs on the contact point. Therefore, the 
deformation of particles is transformed into spring 
deformation under the contact force. Chose two 
similar particles of P, Q in the coarse-grained soil, it is 
shown in FIG. 1. 



FIG 1. THE RIGID CONTACT MODEL OF COARSE-GRAINED 
SOIL 


For particle P, the equilibrium equation can be 
established. 


Where L l (/? ') is the average of contact force of the 
branch vector within n + dii . 


Under the contact force, the strain tensor is £y(x a ), the 
displacement of contact point is u t (x a ) . Ignoring 
rotation of particles, the equation is shown as below 

u,(x a ) = £ lj (x a )L j (x a ,n') ( i,j = x,y ) (4) 

The virtual work of contact force in the unit volume is 

W = ^f i (x a ,n)e ij (x a )L ] (x a ,n') (5) 


Where Factor 2 represents each particle counted twice. 
The virtual work for average stress is 


W = a n £ ii 


From equation (2), (5), considering the symmetry of 
the stress tensor, the equation is obtained as follow. 


1 M 

— Z[ f i (x a ,n)L i G a ,n') 

4V a = 1 


J ' 
-iv 


( 6 ) 


+ f l (x\n)L i (x a ,n')] 
The integral form is 


M 

^ . f Ut 0 n)Lj (n) + fj (d)L, (n)]E(n)dn (7) 


Equation (7) links the macroscopic average stress 
tensor and microscopic first-order tensor together. The 
coarse-grained soil particle size can be classified into 1 
classes, and each class has an average particle size d k 
with Nk particles. Assuming that all the coarse-grained 
soil particle size is d , then 

V s = M 7td 2 , d=- £d k 2 
2 V *=i 

So 


lfr(x a ,n’)(x/ a -x/) 

a = 1 


= Z frcx\n')(xr -x, p ) 


( 2 ) 


m 

7r( 1 + e) 




•\l K [f(9)n j + f J (e)n l }E(e)dd 


( 8 ) 


Where M is the contact number of particle p; n ' is a 

X p -X q I I 

unit branch vector, and n 1 = — j — , I L pq is branch 

\L m \ 1 1 

length. 

Put the summation of all particles into integral form in 
course-grained soils. 

I ft 0 n)Lj (it ')E(H)dn = j fj («)L, (n ')E(n)dQ. (3) 


Where m is the average coordination number; E is 
void ratio. 

Assuming that the contact force is f t (i=n, s; n, s as the 
local coordinate), then the relation between force and 
displacement on the contact point can be expressed as 
in incremental form. 

4/i = DyAU j (9) 
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Where D , is contact stiffness tensor. 

y 

In local coordinates n, s, D iJ can be expressed as 

D j = D n n,n J+ D s S,S J (10) 

Generally, the normal contact stiffness is a function of 
the normal contact force. 

D „=Cjf (11) 


Where C,j3 are function associated with the coarse- 
grained soil characteristic, particle size and surface 
roughness. 

The tangential stiffness is 


D,=C 2 D b ( 1- 


fs 


f„^ncp 




Where C 2 ,r/ are coefficients related to the material 
itself; (j) fl is frictional angle between particles. 


Two-dimensional Granular Material 
Constitutive Model 

Coarse-grained soil particles are randomly stacked, so 
it can be described by fabric density distribution 
function. According to findings of Rothenburg and 
Bathurs, the density distribution fuction of coarse- 
grained soil particle fabric content could be described 
approximately by three fuction in the two-dimensional 
cases. 

Contact normal vector 

E(d) =— (i + ai (cos e-e x )} (12) 

2n 

Normal contact force 

/ B (0) = /o[l + «2 cos 2(0-30] (13) 

Tangential contact force 

/,(*) = -/ 0 «3 sin 2(0- 0 3 ) (14) 

Where the minus is Negative tangential contact force 
to rotate counter clock wise positive; a l ,a 2 ,a 3 are 
coefficients reflecting the degree of anisotropy; 
0 X ,0 2 , 0 3 denote fabric shaft angle, maximum normal 
contact force, average contact angle of maximum 
tangential contact force respectively. 

/ 0 is average tangential contact force. 
fo = lo 

By substitution formula (12-14) into the formula (8) 


<h * = 7 i d l ) ^ t 1 + ^(«i cos 23 + 

(1 + e)zr V*=i J 2 

a 2 cos2 # 2 + a 3 cos 20 3 ) + — (cos 26( cos 2<9 2 (15) 

+ sin 23 sin 23 ,)] 

^22 = 12 [1-T( a i COS 23 

(l + e);rV*=i J 2 

+a 2 cos 23 +a 3 cos 23 ) + ^-^-(cos 23 cos 20 2 (16) 
+ sin 23 sin 26h )] 


<h 2 = n m \ { z d k 1 («1 sin 23 

(l + e)^V<t=i J 

+a 2 sin 20 2 + a 3 sin 20 3 ) 


(17) 


Stress is associated with void ratio, the number of 
contact points and fabric contents. If the relationship 
between contact force and contact displacement is 
linear, the formula is obtained as follow by Hooke's 
law. 


Act = E As 

The tangential stiffness 

D s = 

Where A is a constant, generally A = 0.1 -1.0 . 

The relationship between strain and stress of coarse- 
grained soil can't be derived directly, and the strain 
should be linked to contact force. Coarse-grained soil 
relationships between stress and strain can be derived 
when the local constitutive relation is confirmed. 

The contact stiffness tensor is 

D ij = D n n i n j + D s s i s j (b J = 1> 2 ) (13) 

Local contact constitutive relation of incremental form 
4/; = DyAu j (19) 

Where Af] is incremental contact force on the contact 
point; Aiij is incremental displacement on the contact 
point. 

For the coarse-grained soil, the number of particles is 
homogeneous on the large scale, and its displacement 
is linear distribution without forming shear zone. 

AUj = IjAsy = / (x a )n i As j (20) 

Where Asj is strain increment. 

The relation between contact stress increment and 
strain increment is 
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Af i (x a ,n) = l(x a )(D n n i n J 
+ D s s,Sj)n k As kj 


(*,7 = 1,2) 


( 21 ) 


By substitution formula(21)into formula(20), the 
incremental formula of force and strain is 


relationship between stress and contact force is 
affected by the fabric. The mechanical response of 
coarse-grained soil in the large scale proved that the 
changes of fabric contents play an important role on 
the deformation characteristics of materials. 


a <T/ = Ajki A£ ij ( h j = !’ 2 ) (22) 

Where 


4* = ( w*"/ + WM* (23) 

B ijki = +n J s i n k s l +n i s J n l s k +n j s i n,s k )/ 4 

Density function is 

£(«) = \j2n (24) 

Substitute formula (24) into formula (23), and integral 
is 


Ao_ x 

Ac7 x 

Act 


m 

4(1 + e)^r 


'3 D„+D s 

D n ~ D s 

0 

A Gx 

D n ~ D s 

3 D n +D s 

0 

<1 

0 

0 

D n + D s _ 

1 

> 

1 


(25) 


The rigid contact model of coarse-grained soil is 
established based on the microstructure of coarse- 
grained soil particles. Local constitutive relation of 
coarse-grained soil is acquired by the connection 
between contact force of local stress between particles. 
Based on the local constitutive relation of coarse- 
grained soil, the constitutive relation of two- 
dimensional granular of coarse-grained soil has been 
set up. These studies provide foundations for the 
further researches. 

The application of basic theories of discrete mechanics 
and particle mechanics is a promising method to study 
stress-strain characteristics of coarse-grained soil. The 
two-dimensional constitutive relation has been 
discussed in this paper. For topics related to three- 
dimensional which is quite similar to the two- 
dimensional approaches. And the contact force of 
coarse-grained soil particles and contiguous normal 
distribution need further researches. 


The bulk modulus, shear modulus and Poisson's ratio 
are respectively 

K _ m°n 
2(1 + e)n 

c m(D „ + DJ 
4(1 + e)7i 


3 D n -D, 

It is proved that the bulk modulus and shear modulus 
of Coarse-grained soil is related to pore ratio and 
coordination number. 

Conclusions 

The physical and mechanical properties of coarse- 
grained soil are closely related to particle spatial 
stacking mode, void size and the distribution of void 
space. The relationship between average force and 
contact point is set up by analyzing the force between 
particles and normal direction of contact based on the 
fabric characteristics of coarse-grained soil. Due to 
fabric change during the process of deformation, 
relation between force and strain is nonlinear as same 
as relation between force and contact force, the 
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